AD-A228  749 


FINAL  REPORT 

INVESTIGATION  AND  MODELING  OF 
RADIATION  ABSORPTION  PROCESSES  AND 
OPACITIES  IN  DENSE  PLASMAS 


1990 


FOR  AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 


1401  MCCORMICK  DRIVE,  LANDOVER,  MARYLAND  70785 

(301)  925-9400 

'DISTRIBUTION  statemenTX"  ' 

Approved  for  pnblic  reltOMj 
- DUtribotton  Usllalfrd 


d 


:IT~Y  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION. PAGE 


Form  Approved 
OMB  No.  0704-0188 


:PORT  SECURITY  CLASSIFICATION 
U  CLASSIFIED 


1b.  RESTRICTIVE  MARKINGS 


XURITY  CLASSIFICATION  AUTHORITY 


ECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 
Distribution  Unlimited 


^FORMING  ORGANIZATION  REPORT  NUMg£R(S) 

FA-0184Z 


5.  MONITORING  ORGANIZATION  REPORT  N,U<VttER(S) 

AEOSR-TR-  9  0  1 0  ‘  ^ 


AME  OF  PERFORMING  ORGANIZATION 
FA,  Inc. 


6b.  OFFICE  SYMBOL 
(If  applicable) 


7a.  NAME  OF  MONITORING  ORGANIZATION 

AFuse//uP 


3DRESS  (City,  State,  and  ZIP  Code) 
401  McCormick  Drive 
andover,  MD  20785 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 

!3M<f  mo 

rtolhne.  fiF/3  Pt 


AME  OF  FUNDING/SPONSORING 
RGANIZATION  AIR  FORCE 

0E  OF  SCIENTIFIC  RESEARCH 


8b.  OFFICE  SYM80L 
(If  applicable) 


1 1  p 
r 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F49620-88-C-0055 


10.  SOURCE  OF  FUNDING  NUMBERS 


UILDING  410, 

PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

0LLING  AFB 

ELEMENT  NO. 

NO. 

NO. 

ACCESSION  NO. 

D.C.  20332-6448 

C  Ho  2F 

230  J 

G9 

TLE  (Include  Security  Classification) 

(U) 


INVESTIGATION  AND  MODELING  RADIATION  ABSORPTION  PROCESSES 
AND  OPACITIES  IN  DENSE  PLASMAS 


:RSONAL  AUTHOR(S) 
day  Gupta,  PhD 


TYPE  OF  REPORT 

INAL 


13b.  TIME  COVERED 

from  Q4/01/88toQ2/28/90| 


14.  DATE  OF  REPORT  (Year,  Month,  Day) 
/990 


15.  PAGE  COUNT 

/O  / 


JPPLEMENTARY  NOTATION 


COSATI  CODES 

18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

IELD 

GROUP 

SUB-GROUP 

L3 

BSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 


A  variety  of  atomic  processes  contribute  to  absorption  of  radiation  in  dense  plasmas.  Most  existing  atomic  data 
;  for  low  density,  high  temperature  plasmas.  At  high  densities  and  low  temperatures,  realistic  modeling  to 
:orporate  the  important  additional  effects  in  order  to  generate  atomic  data  is  needed.  This  is  addressed  in  the 
isent  work. 

The  models  and  computer  codes  developed  for  the  project  includes  effects  of  non-linear  screening,  electron 
generacy,  exchange-correlation  and  ion  interactions  self-consistently.  These  were  applied  to  ions  of  specific 
afigurations  in  dense  plasmas  and  represent  improvements  over  'average  atom  models'  often  used  in  dense  plasma 
ysics. 

The  focus  of  this  work  is  mainly  on  the  bound-bound,  bound-free  and  free-free  photoprocesses,  that  contribute  to 
liation  absorption  and  opacity  of  dense,  low  temperature  plasmas.  We  discuss  a  model  to  generate  ionic  distribution 
it  is  computationally  faster  than  rate  equation  method.  We  also  discuss  a  model  to  investigate  the  d.c.  electron 
iduction  in  dense  plasmas,  which  incorporates  effects  of  multiple  scattering  and  improves  over  the  Ziman  type 


>del. 

Continue  on  reverse...  j 

ISTRI0UT.ON  /  AVAILA8ILITY  OF  ABSTRACT 
UNCLASSIFIED/UNLIMITED  53  Z,\MC  AS  HPT. 

□  OTIC  USERS 

I  AM  6  OF  RESPONSIBLE  INDIVIDUAL 

22b.  TELEPHONE  (Include  Area  Code) 
Zo2-  it 

22c.  OFFICE  SYMBOL - 

A/P 

>rm  147 J,  JUN  86 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


Previous  editions  are  obsolete. 


SECURITY  CLASSIFICATION:  Unclassified 


SFA-0184Z 


INVESTIGATION  AND  MODELING  OF  RADIATION 
ABSORPTION  PROCESSES  AND  OPACITIES  IN  DENSE 
PLASMAS 


PREPARED  BY: 
UDAY  GUPTA,  PHD 


SFA,  INC. 

1401  McCormick  Drive 
Landover,  Maryland  20785 


SEPTEMBER,  1990 

CONTRACT  NUMBER  NO:  F49620-88-C-0055 

PREPARED  FOR: 


Accession  For 

liTIS  GRAflbl 
DTIC  TAB 
Unannounced 
Justification. 


* 

□ 


Br - 

Distribution/ 


Availability  Codes 


Avail  and/or 
Special 


,JSL 


AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 
BUILDING  410 
BOLLING  AFB 

WASHINGTON,  D.C.  20332-6448 


1401  MCCORMICK  DRIVE,  LANDOVER,  MARYLAND  20785 

(301)  925-9400 


TABLE  OF  CONTENTS 


Page 


INTRODUCTION 

I  ATOMIC  MODEL  FOR  IONS  OF  SPECIFIC  CONFIGURATIONS  IN 
PLASMAS 

H  BOUND -BOUND  PHOTOPROCESS 
IH  BOUND-FREE  PHOTOPROCESS 

IV  FREE-FREE  TRANSITIONS 

V  MODEL  FOR  DISTRIBUTION  OF  IONIC  SPECIES  IN  PLASMAS 

VI  ELECTRICAL  CONDUCTIVITY  IN  HOT  DENSE  PLASMAS 
RESULTS 

CONCLUSIONS 

ACKNOWLEDGEMENTS  <  .  ’  ' 

V  " 

'  c  •  ' 

REFERENCES  '"V 

APPENDIX  I 

FIGURE  1  '  ..  • 

/ 

'  7 

TABLES 

SAMPLE  COMPUTER  OUTPUT 

LISTING  OF  SAMPLE  COMPUTER  CODES:  FFNEW.FOR;  DIFF.FOR; 
SIMP. FOR;  WAVEF.FOR. 


1 

4 

6 

8 

8 

10 

15 

19 

22 

22 

23 

24 


INTRODUCTION 


Modeling  of  radiation  absorption  processes  in  plasmas  require 
investigation  and  computation  of  a  large  number  of  atomic  processes 
contributing  to  such  absorption.  Bound-bound  phototransition  (or 
photoexcitation),  bound-free  (or  photoionization)  and  free-free  (or  inverse 
Bremsstrahlung)  photoprocess,  as  well  as  Thomson  and  Compton  scattering  of 
radiation  and  electron  conduction  of  low  energy  radiation  are  some  of  these 
processes.  This  report  deals  with  the  investigation,  modeling  and  computation 
of  these  atomic  processes  in  plasmas. 

Depending  on  the  plasma  density  and  temperature,  several  different 
effects  need  to  be  systematically  investigated  and  modelled.  The  bulk  of 
theoretical  and  experimental  atomic  data  available  in  the  literature  is  limited  to 
isolated  neutral  atoms.  Depending  on  experimental  conditions,  plasma  ions 
can  be  ionized  to  a  different  degree  of  ionization.  The  atomic  data  for  highly 
charged  ions,  as  occur  in  high  temperature  Takamak  plasmas,  astrophysical 
plasmas1  or  X-ray  laser  plasmas  are  quantitatively  very  different  from  those  for 
neutral  atoms.  Also,  at  higher  plasma  densities  (and  for  comparatively  lower 
temperatures),  plasma  screening  effect  causes  ionic  energy  levels  to  shift 
significantly  (particularly  the  upper  levels)  and  modifies  the  effective  ionic 
potentials  so  that  the  various  atomic  processes  can  be  significantly  altered. 
Atomic  data  for  isolated  atoms  or  ions  are  no  longer  valid  for  these  dense 
plasma  conditions.  Realistic  modeling  of  atomic  processes  are  therefore 
required  for  dense  plasmas.  In  addition,  broadening  of  ionic  spectral  lines  also 
occur  at  higher  densities  (such  as  in  ICF  plasmas)  which  influence  radiation 
absorption  processes.  Line  broadening  is  temperature  sensitive  as  well. 

At  extremely  high  plasma  densities,  such  as  in  superdense  ICF  pellets  or 
in  astrophysical  plasmas  (interior  of  White  Dwarfs,  for  example),  processes 
such  as  Compton  scattering  of  radiation,  as  well  as  electron  conduction  of 
radiation  becomes  important  for  overall  absorption  processes  and  opacities. 

For  plasmas  in  those  conditions,  these  processes  therefore  need  to  be 
incorporated  in  the  modeling. 
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The  modeling  and  computation  of  radiation  absorption  coefficients  and 
plasma  opacities  are  required  inputs  to  overall  radiation  hydrodynamic 
simulation  of  plasmas  in  a  variety  of  laboratory  and  astrophysical  conditions. 
We  are  concerned  with  low  to  fairly  high-density  plasmaS  as  occur  in 
laboratory  discharge  tube  plasmas,  Z-pinch  plasmas,  Takomak  plasmas.  X-ray 
laser  plasmas,  wire-imploded  plasmas,  ICF  and  laser  imploded  plasmas,  and  for 
some  astrophysical  plasmas,  such  as  outer  layers  of  White  Dwarfs.  The  atomic 
data  for  radiation  absorption  are  essential  ingredients  for  radiation  transport 
simulations,  as  well  as  overall  radiative  (and  spectroscopic)  studies  of  plasmas 
in  above  mentioned  experimental  conditions. 

The  effect  of  plasma  density  and  temperature  on  atomic  processes  has 
been  treated  in  the  past  by  'average  atom  models'  (AAM's)  4  (The  Thomas- 
Fermi-Dirac  model13d4  js  one  Gf  the  -widely  used  AAM's).  However,  there  are 
well-known  serious  shortcomings  of  the  AAM's.  These  models  assign 
fractional  occupation  of  bound  electronic  states  and  thus  lead  to  fictitious 
atomic  structure.  The  AAM's  do  not,  therefore,  generate  accurate  data. 

For  accurate  calculations  of  atomic  processes  for  ions  in  dense  plasmas, 
the  correct  treatment  of  plasma  density  and  temperature  effects  on  the  atomic 
structure  of  multielectron  ions  of  specific  configuration  is  of  central 
importance.  The  specific  ionic  configurations  with  integral  occupation  of 
bound  electronic  states  need  to  be  considered  in  order  to  obtain  realistic  atomic 
structure.  Screening  by  plasma  electrons  has  to  be  treated  in  a  fully  non-linear 
fashion  (since  Debye-Huckel  screening  is  known  to  be  grossly  inaccurate  for 
dense  plasmas).  Screening  by  both  bound  and  free  electrons  has  to  be 
considered  in  a  self-consistent  way.  Degeneracy  effect  and  many-electron  effect 
such  as  exchange-correlation  interaction  need  to  be  accounted  for  in  the  model. 
Depending  on  plasma  conditions,  the  effect  of  strong-coupling  (i.e.  when 
Coulomb  interaction  dominates  over  random  thermal  motion)  needs  to  be 
taken  into  account. 

The  main  objective  of  the  present  work  is  to  incorporate  these  effects 
systematically  in  the  models  and  the  relevant  computer  codes  in  order  to  have 
the  capability  of  generating  realistic  atomic  data  for  various  processes 
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contributing  to  photon  absorption.  The  bound-free  (photoionization),  the  free- 
free  (inverse  Bremstrahlung)  and  the  bound-bound  (photoexcitation)  atomic 
processes  arc  of  particular  interest  in  this  context.  The  improved  models  and 
the  various  computer  codes  would  be  useful  tools  for  such  calculations  later  on 
particularly  for  moderately  high  electron  density  range  (~1023-1024  e/cm3)  and 

relatively  lower  temperature  range  (~10-200  eV . )  where  'average  atom'  type 

models  are  known  to  be  inadequate. 


I.  ATOMIC  MODEL  FOR  IONS  OF  SPECIFIC  CONFIGURATIONS  IN 

PLASMAS 

For  a  given  ion  immersed  in  a  plasma  of  mean  electron  density  pe  and 
mean  ion  density  p:  at  temperature  T,  we  have  a  coupled  system  of  equations 
for  electrons  and  ions  to  solve.  For  electrons,  we  have, 

-4.v*+  VtCr)]4>.<r)  =  HD 


pe(r)=  Z  fbHHCr)!2- 

fc> 


and  for  ions 


_  -v^/fegT 

ft  e 


i-pn 
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with 


These  coupled  equations  are  solved  self-consistently.  In  equation  1(1)- 
1(5),  pc(r)  and  p;  (r)  are  the  electron  and  ion  distribution,  Z  is  the  mean  ionic 
charge  (Z=  Z  -  nb)  the  number  of  bound  electrons.  R  is  the  radius 
(approximately  the  ion  sphere  radius)  over  which  charge  neutrality  is  achieved. 
/' s  are  the  occupation  factors.  For  ion  of  a  particular  configuration, /h's  are 
integral  occupation  factors  for  the  bound  states  of  the  ion.  In  contrast;  lor 
’average  atom  model',  the  f's  are  taken  to  be  Fermi  distribution  functions  with 
fractional  occupation,  qb's  are  the  bound  and  <hkl’s  are  the  continuum 
wavefunctions  of  energy  k2/2  and  angular  momentum  1.  The  exchange- 
correlation  potential  is  denoted  by  Vxc.  It  depends  on  the  plasma  density  and 
the  temperature  and  is  calculated  ^/3  separately  by  summing  up  the  exchange 
and  "ring"  diagram  sums  and  then  used  as  input  to  subsequent  calculations  in 
eq.  I(l)-I(5). 

Let  us  point  out  that  the  effective  ionic  charge  Z  (=  Z  -  nb,  nh  is  the 
number  of  bound  electrons)  is  an  essential  ingredient  to  the  calculation.  In  our 
model,  an  'average  atom'  model  -type  calculation  is  done  first  to  generate  Z . 

For  this;  equations  I-(l)  through  I-(5)  are  solved  numerically  to  self-consistency 
with  fb's  and  f(k)’s  in  eq.  I-(3)  taken  as  Fermi  functions  (i.e.  fractional 
occupation  numbers). 

Once  Z  is  generated  from  self-consistent  'average  atom'  type  model  as 
outlined  above,  it  is  then  used  as  input  to  subsequent  iterative  calculation. 
Thus,  eq.  I-(l)  through  I-(5)  are  solved  numerically  to  self-consistency  with  f,,'s 
taken  to  be  integral  occupation  factors  corresponding  to  the  ion  of  specific 
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configurations,  (for  example,  for  Ne-like  ground  state:  IS2  2S2  2p6;  fb 's  are  2 
for  IS-state;  2  for  2S-state  and  6  for  2p-state).  This  is  the  important  difference 
between  our  model  for  ions  of  specific  configuration  and  the  'average  atom" 
type  model4. 

Once  the  self-consistent  potentials  for  ions  in  plasmas  are  obtained,  the 
bound-free,  bound-bound  and  free-free  cross  sections  can  be  computed  by 
generating  the  relevant  numerical  bound  and  continuum  wavefunctions  first 
and  then  computing  the  relevant  matrix  elements.  These  photoprocesses  will 
be  outlined  in  the  next  section. 

n  BOUND-BOUND  PHOTOPROCESS 

For  absorption  of  a  photon  of  energy  fia>  in  a  spectrum  line  so  that  an 
electron  makes  a  transition  from  one  ionic  state  xa  to  another 

M  — *•  %b  (II  I) 

the  photon  frequency  is  close  to  the  line-center  frequency  co0: 

k«0  =  Eh-Ea  <"-2> 

The  contribution  to  the  monochromatic  opacity  Kwdue  to  (II-l)  is: 

Kj^b)  =  NUf Cm)  4>m  «i-3> 

me 

where  Na  is  the  number  density  of  ions  Xa,  is  a  profile  factor  normalized  to: 
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(II-4) 


and  f  (b,a)  is  the  oscillator  strength  given  by9: 


f  (b,a)  = 


where  S  (b,a)  is  the  line  strength: 


S  Q.  a.) 


S(b,a)=|<bllDl|a.>|Z 


(II-5) 


(H-6) 


with  the  operator  D  is  dipole  operator. 

If  the  bound-state  wave  functions  (numerically  generated  from  the  self- 
consistent  potential)  are  represented  as: 

\  (?)  =  (r)  ^ 

where  Yim(?)  is  the  spherical  harmonic.  Then, 

S  ( t,  i')  =  Z  KJ™  I  ?  1 

m  w' 

=  Z 

‘Tv'n,/<V  ^  (II-8) 

=  l>  [J  ft  Pi'  «*  «*r] 

Where  1,T  are  the  angular  momentum  values;  m,m'  are  the  Z  components  and 
1>  is  the  greater  of  the  1  values  considered. 
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m  BOUND-FREE  PHOTOPROCESSES  (PHOTOIONIZATION) 


For  photoionization; 

we  replace  \| /IW  by  yVm 

,  defined  by: 

Hiv  = 

ti  v.-  w  i 

(III-l) 

with 

Fkt  « 

^  — - —  Sin  ( 

co  V^'/ 

K  +  ■  •  ■') 

(ffl-2) 

for  normalization  condition  for  \j/r  ,is: 

=  1 

(III-3) 

and  /  ^ 

<%..  I 

*«■>  - 

S(e-e') 

The  oscillator  strength  for  photoionization  is  given  by: 


df  A  1 0 

3n  (21 +i) 


Iff 

V 


k't' 


(III-4) 


IV  FREE-FREE  TRANSITIONS 


To  obtain  the  oscillator  strength  df/dE  for  photoabsorption  of  one 
photon  co  in  the  field  of  an  ion  in  a  plasma  of  electron  density  Ne  of  mono- 
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energetic  electrons,  we  replace  in  the  expression  for  bound-bound  by 
and  VimbyVimk  'vhich  represents  the  lm  component  of  the  total  function 

describing  the  flux  of  eiectrons  having  momentum  k.  For  isotropic 
distribution  of  velocities,  such  function  has  the  asymptotic  form 


i  Ef*r 

e  (if  rV-*o)  (IV'1} 


(and  an  analogical  expression  for  Coulomb  Field). 

It  can  be  shown  that,  the  oscillator  strength  for  the  free-free  process  is  given  by: 
(see  Appendix  I  for  a  proof) 
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V 


MODEL  FOR  DISTRIBUTION  OF  IONIC  SPECIES  IN  PLASMAS 


For  calculation  of  overall  radiation  absorption  and  opacities,  it  is 
necessary  to  know  the  distribution  (i.e.  number  densities)  of  ions  of  different 
types  that  occur  in  plasma  at  a  particular  density  and  temperature.  It  is  possible 
to  obtain  these  distributions  from  detailed  rate  equation  calculations.  The 
computation  involved  in  rate  equation  solution  is  time-consuming,  complex 
and  tedious.  Here  we  discuss  an  alternative  approach.  This  model  is  based  on 
the  density  functional  method,  following  Perrot.8  This  outlined  below  to 
point  out  specific  aspects  of  the  calculation  involved. 

Opacities  of  plasmas  have  often  been  calculated  using  'average  atom’ 
model  (AAM)  for  ionic  species.  The  AAM,  however,  deals  with  fictitious 
’average  atom’  with  fractional  occupation  of  electronic  bound  states.  In  a  real 
plasma,  a  number  of  different  ionic  species  are  usually  present  in  various 
abundance. 

Let  us  first  consider  the  total  energy  of  an  ionic  species  defined  by  the 
configuration  of  the  bound  spectrum  and  a  screening  charge  of  free  electrons  in 
LTE.  The  fixed  integer  occupation  numbers  of  bound  electrons  are  n^s  while 
the  fractional  occupation  numbers  of  free  electrons  are  denoted  by  fk’s.  The 
total  energy  considered  here  is  that  of  embedding  the  atom  in  the  plasma: 


AE  ( - 8=3  E  (.  ...  Tlj,  ,••••)-  I  [  p] 


(V-I) 


where  E(p)  is  the  unperturbed  electron  gas  energy.  One  can  show  that8 

ArjL  •  'l 


oo 


AE 

*  O 

<  f  i™*  c  ft + ft  -  P>>  - 1  <(&-*■  ft  -  PH 0  Cfi-A- 

C  ( fb)  -  <  fb  LVx°,  ( fO  +  wxct  (p, ,  f  )]>  + 
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CCfti  -  C  c  F)  -  <(V vx7  ( ft)>  -  vxl(r)> 


-1% 


(V-2) 


^  R 


In  (V-2);  the  notations 

4f*> 

f-9 

and 


J f  (r) 

J  fCr ')  $(r -?')*?' 


(convolution  product) 


have  been  used. 

In  (V-2),  pb  is  the  bound  and  pj  the  free  electron  density,  rij  the  phase  of 
the  continuum  electronic  states,  p  the  mean  electron  density  of  the  plasma;  E[ 
is  the  exchange-correlation  energy,  and  Wxcb  is  the  bound-free  exchange- 
correlation  potential. 

Next,  consider  the  transition  of  an  electron  from  initial  state  i  to  final 
state  j  in  the  bound  spectrum.  The  change  in  the  total  energy  of  the  species 
associated  with  the  transition  is: 


SEZ„  C  ‘ 


j)  =  AE.  (.  .. 

—  AE  (- 


In  above,  Z*  refers  to  the  ionization  degree  of  the  species.  Energies  in  eq. 
(V-3)  are  defined  by  cq.  (V-l).  The  ’transition-state’  theory  of  Slater  is  utilized 
in  the  calculation  here.  It  improves  on  the  AAM  which  has  been  used  often  to 
obtain  the  average  degree  of  ionization  in  a  plasma. 

In  a  plasma  of  different  ionic  species,  let  us  consider  the  most  abundant 
species  with  a  given  ionization  degree  Z*.  The  energy  of  this  one,  which  has 
the  configuration  <...ni()...>  is  chosen  as  a  reference.  The  total  energy  of  another 
species  with  the  same  Z*  but  a  different  configuration  (...iij...)  can  be  shown  to 
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be: 

A£-g>  (*  •••  nu:)  tz*  +  Q  (An'2-') 


(V-4) 


where  Xeiz*  is  the  modified  effective  energies  of  the  referenced  species, 
provided  that 


2.  An;  =2Xnt  -Mio)  =  0 

bound  bound 

The  system  can  now  be  interpreted  as  a  collection  of  independent 
pseudoparticles,  distributed  among  states  of  energy  eiz*.  The  average 
occupation  numbers  for  the  species  will  be  given  by  FD  statistics. 


>  +  **pL?C 


(V-6) 


pz*  is  our  effective  chemical  potential  for  the  bound  pseudoparticles,  to  be 
determined  from 


22^1  =  Z  -  Z*  .  (V-7) 

Equations  (V-6)  and  (V-7)  thus  fully  define  the  average  species  for  the 
ionization  degree  Z*.  Finally,  the  probability  of  a  Z*  species  is  given  by 


T|*  =  A  exp  (-  p  Fz  *)  (v-8) 

where  Fz*  is  the  effective  free  energy  which  can  be  obtained  from  the  AE^, 
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(...n;,...)  discussed  above  and  the  entropy  of  the  bound  pseudo-particles. 


sz *.  =  -  kb  £.[?  ‘ 1,1  (v'9> 

In  (V-8),  A  is  a  normalization  constant  such  that  Epz*  =  1,  the  sum  running 
over  all  integer  ionization  states, 

and 

Fz,  =  ABZ,  (•  •  •  .nt,  -)  _TSZ*  -  T*S  [ft] 

(V-10) 


is  the  effective  free-energy  determining  the  probability  via  eq.  (V-8). 

As  an  example  of  the  application  of  this  model^,  consider  probabilities  of 
various  ionic  species  in  the  Aluminum  plasma  at  temperature  100  eV  and 
electron  density  of  4.8  x  1023  cm'3.  The  average  atom  model  gives  the  average 
ionization  degree  for  the  plasma  in  the  above  condition  as  7.9.  This  indicates 
that  the  effective  ion  charge  Z*  =  8  (i.e.  ions  with  5  bound  electrons  or  Boron¬ 
like  ions)  will  have  the  highest  probability  of  occurrence.  This  is  confirmed 
from  the  calculation  using  the  above  model8 :  Pz*  is  highest  for  Z*  =  8  and 
drops  off  rapidly  for  lower  and  higher  Z*  ionic  species.  The  relative 
abundances  of  ions  of  different  species  in  the  Al-plasma  at  the  above  density 
and  temperature  are  thus  obtained. 

It  is  seen8  that  for  the  above  plasma  conditions  B-like  A1  is  most 
abundant,  followed  by  Be-like,  C-like,  O-like,  Li-like  ions  etc.  in  decreasing 
order  of  relative  abundances.  Thus  the  model  discussed  above  provides  an 
useful  tool  to  generate  the  number  densities  (or  abundances)  of  different  ionic 
species  in  a  plasma.  These  data  are  necessary  for  subsequent  computation  of 
radiation  absorption  in  a  plasma. 

Note  that  the  distribution  of  ionic  species  in  plasmas  are  usually 
obtained  by  solution  of  rate  equations.  The  rate  equation  method,  however,  is 


-  13  - 


very  complex,  tedious  and  time-consuming.  In  view  of  the  advantages  of  the 
model  discussed  above,  we  are  in  the  process  of  building  a  computer  code  in 
order  to  generate  the  ionic  distribution  in  different  plasma  conditions. 
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VI  ELECTRICAL  CONDUCTIVITY  IN  HOT  DENSE  PLASMAS 


Electron  conduction  is  important  in  absorption  of  radiation  from  lasers 
because  low  frequency  inverse  bremsstrahlung  is  essentially  the  same  as  Ac 
joule  heating.  Emission  of  visible  and  ultraviolet  light  (and  generation  of 
magnetic  field)  also  involve  electron  current  flow.  Electron  heat  conduction  is 
an  important  limitation  on  thermonuclear  burn  and  is  particularly  important 
in  laser-target  ablation.  The  energy  transported  by  electron  conduction,  in  turn, 
determines  conductive  opacity. 

In  this  section,  we  discuss  a  model  to  compute  electrical  conductivity  (or 
resistivity)  of  dense  plasmas.  Various  methods,  such  as  Kubo  approach  or 
Boltzmann  transport  eqn.  method  has  been  used  by  different  groups  to  obtain 
electrical  conductivity'7  (or  resistivity)  in  plasmas.  A  successful  approach  is  via 
Ziman  formula^  which  is  essentially  a  variational  solution  to  Boltzmann  eqn. 
Following  Perrot  and  Dharma-Wardana^,  the  method  we  will  discuss  here  is  a 
generalization  of  Ziman  formula  using  density  functional  method.  This 
extends  the  applicability  from  the  usual  weak-isolated-scatterer  limit  to  the 
regime  of  strong  multiple  scatterers,  as  is  important  for  high  density  plasmas. 

a)  REGIME  OF  ISOLATED  SCATTERERS 


The  most  familiar  form  of  Ziman  formula  applies  to  plasmas  and  liquid 
metals  in  regime  of  weak  and  strong  but  isolated  scatterers.  The  resistivity  R  is 
given  by. 


J _ 

Znz2*?-? 


OQ  2K 


(VI-1 ) 


where  q2  =  2k2  (1-Cos  0).  Here  q  is  the  momentum  transferred  from  the 
incident  electron  with  energy  e  =  k2 .  The  derivative  of  the  Fermi  distribution 
function  for  electrons  at  a  mean  density  n  =  Z  p  is  denoted  by  f(e).  The  mean 
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ion  density  is  p  .  The  ion  distribution  is  specified  by  the  structure  factor  S(q). 
The  differential  scattering  cross-section  c(q)  depends  on  the  incident  electron 

momentum  k  and  the  transferred  momentum  q.  For  strong  scatterers  o(q)  has 
to  be  obtained  from  the  phase  shifts  of  Vci. 

In  the  weak  isolated  scatterer  limit,  the  Born  approximation  may  be  used 
to  obtain  c(q): 

o(q)  =  IVei  (q  )  /  4rc  e(q)|2  (VI-2) 

In  above  Vci  (q)  is  the  weak  electron-ion  interaction,  and  e  (q )  is  the  exact 
interacting  uniform  electron-gas  dielectric  function  at  a  density  n  and 
temperature  T.  The  density  functional  method  outlined  previously  is  capable 
of  generating  e(q,n,T.)  also. 

The  ion-ion  structure  factor  S(q)  (VI-1)  is  determined  by  the  ion  density 
p,  temperature  T,  and  the  ion-ion  pair  interaction  V;j  (q  ).  In  the  DFM  model  of 
calculation  for  single  ions  immersed  in  an  electron  gas,  the  Kohn-Sham 

equations  are  solved  numerically  for  continuum  electron  states  <hkl  (r )  which 
incorporates  the  appropriate  phase  shifts  5,  (k).  At  self-consistency,  the  phase 
shifts  satisfy  the  finite  temperature  version  of  the  Friedel  sum  rule,  thus 
assuring  the  correct  construction  of  the  continuum  electron  density  of  states  as 
modified  by  the  scattering  potential.  The  scattering  cross-section  is  given  in 
terms  of  the  phase  shifts  by: 

<r(t.e)  =  Vco  ^ 

Y~k  Z(2U  0  (e  ^  — l)  6Kti)fvi-3> 

The  scattering  angle  0  and  the  momentum  transferred 
q  =  k'  -  k  is  related  by 

cfr  =  2.K7-  (l-C-o-se)  (VI-4) 
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2 

For  elastic  scattering  k2  =  k'  =  e. 


b)  THE  REGIME  OF  STRONG  MULTIPLE  SCATTERERS 

The  ion-ion  interactions  in  this  regime  cannot  be  treated  as  a  sum  of 
independent  scattering  centers,  either  for  weak  or  strong  scatterers.  Thus,  the 
electron  is  scattered  from  a  given  quasi-static  ion  distribution  consisting  of 
many  centers.  For  a  definite  configuration  C  of  ions,  the  resistivity  can  be 
written  as  : 


OO 

-  Ar,  f'«>  Fn3  tkA> 

3h  e  va  J  J  ^ 

a 


(VI-5) 


In  the  DFM  method,  the  T-matrix  for  the  given  a  configuration  C  is 
approximated  by  Tkk,  where  Tkk-(e)  is  calculated  from  the  average-ion 
distribution  p  (r).  The  total  T-matrix  of  the  average  distribution  Tkk^e)  is 
obtained  from  the  phase  shifts  Aj  (k )of  the  Kohn-Sham  equation  solved  in  the 
external  field  of  the  scattering  ion  and  its  associated  distribution  p  (r).  In  the 
isolated  scatterer  limit. 


(£)  —  f  <1  *-*’0  *0*- *  '  0  (VI-6) 

In  general,  the  modified  fermi  function  f  in  eqn.  (VI-5)  needs  to  be 
calculated  from  one-particle  Green's  function  with  self-energy  correction. 
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Perrot  et  al^  have  shown  that  eqn.  (VI-5)  can  be  simplified  to 


Ja _ 

3tt  e  n 


Ut  fee-) 


ZK 


H  f3  t(<0 


o 


o 


(VI-7) 


for  strong  multiple  scatterers. 

It  is  clear  from  above  discussion  of  the  model  that  the  structure  factor 

S(k)  and  the  T-matrix  of  the  average-ion  distribution  Tkk-  (e)  are  crucial 
ingredients  in  the  calculation  of  resistivity  of  dense  plasmas  inclusive  of  strong 
multiple  scattering  effect.  In  our  self-consistent  calculation,  we  can  generate 
the  structure  factor  by  taking  fourier  transform  of  pair  correlation  function 
which,  in  turn,  is  generated  from  charge  distributions  obtained  self- 
consistently.  The  phase  shifts  Aj  (k)  in  the  field  of  the  scattering  ion  and  its 
associated  distribution  p  (r)  can  also  be  generated  by  our  atomic  codes.  We  are, 
therefore,  in  a  position  to  utilize  the  model  discussed  above  to  calculate 
electron  conductivity  inclusive  of  multiple  scattering  effects.  Let  us  note  that 
the  energy  transported  by  electron  conduction,  in  turn,  determines  conductive 
opacity.  However,  conductive  opacity  dominates  over  radiative  opacity  only  at 
super  high  densities  (i.e.  for  extremely  degenerate  plasmas).  For  most 
laboratory  plasmas,  conductive  opacity  is  small  compared  to  radiative  opacity. 
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RESULTS 


The  calculations  for  the  atomic  processed  contributing  to  absorption  of 
radiation  involve  development  of  several  computer  codes.  A  large  part  of  the 
work  for  the  project  involved  development,  modification,  testing  and 
benchmarking  a  set  of  computer  codes  for  different  plasma  conditions. 

For  the  project,  the  following  set  of  codes  were  either  set  up  or  modified 
according  to  the  models  discussed  in  the  preceding  sections: 

1.  An  'average  atom’  code  for  dense  plasmas  according  to  tire  prescription 
of  finite  temperature  density  functional  method^  (This  is  verv  similar  to  the 
model  used  by  Liberman^). 

2.  An  atomic  structure  code  for  isolated  multielectron  ions  of  specific 
configurations. 

3.  A  self-consistent  atomic  structure  code  for  multielectron  ions  erf  spent ic 
configuration  immersed  in  a  plasma.  This  includes  plasma  screening  and 
other  effects  discussed  in  Section  I. 

4.  A  code  to  compute  bound-free  photoionization  cross-sections  for  isolated 
ions.  (Also  a  relativistic  photoionization  code  applicable  for  high  Z-elements). 

5.  A  code  to  compute  photoionization  cross-sections  of  ions  of  specific 
configuration  inclusive  of  dense  plasma  effects. 

6.  A  self-consistent  code  to  generate  free-free  photoprocesses  for  ions  of 
given  configuration  inclusive  of  plasma  effects.  Full  numerical  continuum 
wavefunctions  in  self-consistent  potentials  are  used  in  calculating  matrix 
elements  for  free-free  transitions  in  this  code. 

7.  A  code  to  compute  bound-bound  photo-excitation  cross-sections  for  ions 
of  specific  configurations  inclusive  of  plasma  effects. 

For  a  given  plasma  condition  (i.e.  at  a  given  density  rp,  and  temperature 
T,  first  a  self-consistent  calculation  using  'average  atom'  code  is  done  to  obtain 
the  effective  or  'average'  ionic  charge  Z.  This  is  then  used  in  subsequent 
calculations.  For  an  ion  of  specific  configuration  immersed  in  the  plasma  of 
the  above  density  ne  and  temperature  T,  the  self-consistent  code  described  in 
Sec.  I  is  used  to  generate  energy  levels,  wavefunctions,  self-consistent  ionic 
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potentials  and  other  quantities.  Then  calculations  for  bound-bound 
photoexcitation,  bound-free  photoionization  and  free-free  photoprocesses  are 
done  utilizing  the  numerical  wavefunctions  etc.  The  free-free  cross-sections 
for  an  ion  are  calculated  for  various  angular  momentum  states  (1-values)  for  a 
given  photon  energy.  The  allowed  transitions  are,  of  course,  between 
continuum  states  of  angular  momentum  1  and  1±1).  All  the  contributions 
from  different  1  values  are  summed  up  to  obtain  the  resultant  oscillator 
strength  of  df/dE. 

For  purpose  of  illustration,  we  now  discuss  some  representative  results. 
Results  are  tabulated  for  self-consistent  calculations  of  aluminum  plasma  at 
electronic  density  n  =  2.2xl()23  cnv3  and  temperature  T  =  50  eV.  Computed  data 
for  bound-free  oscillator  strengths  (Ryd'1)  for  photon  energies  in  the  range  0.01 
-  100  Ryd.  are  tabulated  in  Tables  1-3  for  Ne-like  ion  (IS2  2S2  2p6)  in  the  plasma 
in  above  condition.  Table  1  is  for  photoionization  from  lS-state,  Table  2  for  25- 
state  and  Table  3  for  2p-state.  The  typical  variation  of  bound-free  oscillator 
strength  is  seen  to  be  a  rapid  rise  from  threshold,  going  through  a  peak,  then  a 
sharp  decrease  followed  by  a  gradual  decrease  with  increasing  photon  energy. 
This  is  seen  clearly  from  Fig.  1.  Results  for  other  types  of  ions:  IS2  2S2  2p5 
(Tables  4-6);  IS2  2S2  2p4  (Tables  7-9);  IS2  2S2  2p6  3S  (Tables  10-12)  in  Aluminum 
plasma  at  electron  density  nc  =  2.2xl023  cnv3  and  T  =  50  eV  over  the  range  of 
photon  energy  0.01  -  100  eV  are  also  tabulated. 

The  oscillator  strengths  for  free-free  photoprocesses  for  Ne-like  A1  in  an 
Aluminum  plasma  nc  =  2.2xl023  cm  3,  T  =  50  eV  for  a  set  of  photon  energies  are 
given  in  Tables  13-20.  Full  numerical  continuum  wavefunctions  in  self- 
consistent  potentials  are  used  in  evaluating  the  allowed  1  — >(1±1)  transitions 
and  eqn.  (IV-2)  is  used  in  calculating  the  matrix  elements  for  oscillator 
strengths  for  a  range  of  incident  electron  energies.  This  method  of  calculating 
free-free  cross-sections  clearly  is  an  improvement  over  semi-classical 
approximations  such  as  Born-El wert  method12.  For  purpose  of  comparison, 
results  of  free  free  calculations  using  unscreened  Coulomb  potential  for  Ne-like 
A1  (instead  of  self-consistent  screened  potential)  are  also  given.  The  computer 
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output  for  free-free  cross  sections  display  the  relevant  matrix  elements  for 
different  1  values.  For  each  photon  energy,  oscillator  strengths  are  shown  in 
the  computer  printout  for  a  set  of  incident  electron  energies  Ejn's.  The  Ein's 
considered  are  1.0,  2.0,  5.0,  10.0,  20.0,  30.0,  50.0  at  units.  The  large  difference  in 
oscillator  strengths  for  the  same  photon  energies  for  self-consistent  results  with 
those  of  unscreened  Coulomb  potential  show  the  importance  of  including 
plasma  effects  appropriately  in  the  calculations. 

As  an  example  of  a  typical  bound-free  calculation,  a  sample  computer 
output  of  photoionization  from  lS-state  of  Ne-like  A1  in  A1  plasma  is  attached. 
The  self-consistent  potential  is  supplied  as  input.  The  output  shows  some  of 
the  radial  mesh  points  and  the  corresponding  potential,  the  iterative  process 
(by  Numerov  method)  of  solving  the  energy  eigenvalue  of  the  lS-state.  For 
each  photon  energy,  the  corresponding  bound-free  matrix  element  and  the 
cross-section  is  computed  and  tabulated.  The  data  corresponds  to  those  of 
Table  1 . 

As  an  example  of  the  codes  used  for  the  calculations,  the  listing  of 
FFNEW.FOR  for  generating  free-free  oscillator  strengths  for  (1+1) 
transitions  is  also  attached.  Again,  self-consistent  potential  for  an  ion  in  a 
plasma  is  fed  in  as  input  read  by  subroutine  VREAD.  OMEGA  is  a  typical 
photon  energy  in  atomic  units,  FFOSC  is  the  oscillator  strength  summed  over 
all  continuum  states.  Numerical  derivative  of  the  potential  is  used  in  the 
calculation  of  matrix  elements.  The  associated  subroutines  DIFF.FOR  and 
SIMP.FOR  are  also  attached.  DIFF  generates  the  numerical  continuum 
wavefunctions.  Minor  changes  are  needed  in  FFNEW.FOR  for  calculating 
l->  (1-1)  transitions  for  different  photon  energy  OMEGA.  The  incident  electron 
energy  EIN  is  read  in;  EOUT  is  the  outgoing  electron  energy  (EOUT  =  EIN  + 
OMEGA). 

For  overall  calculations  of  opacities,  other  effects  such  as  line  broadening 
and  scattering  of  photons  need  to  be  properly  included  in  the  calculation.  We 
are  in  the  process  of  building  up  the  necessary  computer  codes  to  combine  these 
effects  in  the  overall  computation  of  radiative  opacity  of  dense  plasmas. 
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CONCLUSIONS 


A  number  of  self-consistent  models  and  computer  codes  have  been 
developed  which  are  very  useful  for  investigating  the  large  number  of  atomic 
processes  contributing  to  radiation  absorption  in  dense  plasmas.  The 
important  effects  arising  in  dense  plasmas  --  non  linear  screening,  electron 
exchange-correlation  and  degeneracy,  strong  ion  coupling  etc.  are  systematically 
included  in  these  models.  These  are  applicable  to  ions  of  specific 
configurations  in  dense  plasmas  and  thus  represent  improvements  over 
average  atom  models.  These  methods  are  particularly  useful  in  moderately 
high  densities  (~1023  -  1024  e/cm3)  and  relatively  lower  temperature  range 
(—10  eV-200  eV)  where  average  atom  models  are  known  to  be  inadequate. 

A  large  number  of  atomic  processes  contribute  to  radiation  absorption  in 
plasmas.  Most  of  the  work  for  this  project  dealt  with  bound-bound,  bound-free 
and  free-free  photoprocesses  for  ions  in  dense  plasmas.  The  self-consistent 
calculations  are  involved  and  time-consuming  (particularly  since  the  VAX- 
8350  is  the  only  computer  we  could  use  for  this  project  and  we  did  not  have 
access  to  Cray,  for  example).  For  overall  calculation  of  opacities  in  wide  range 
of  plasma  conditions,  vast  atomic  data  bases  need  to  be  generated  involving 
enormous  amounts  of  computation.  We  hope  to  address  those  tasks  in  future 
years. 
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APPENDIX  I 


FREE-FREE  TRANSITIONS 


For  deriving  eqn.  (IV-2)  we  write  the  wave  function  as: 

?=  M  z(n+<)iV^)  **  ^ 
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-IS 
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Therefore 
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Total  oscillator  strength  per  ion  is  equal  to  the  sum  of  all  possible 
contributions: 


Vi'mV  ^oes  not  depend  on  direction  k'  of  the  scattered  electron,  the  averaging 
over  angles  gives  just  the  form  of  used  here,  co  is  the  photon  frequency. 
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Since 


|  jXl  Fk'l' 

we  have 
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This  proves  eqn.  (IV-2). 
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FIGURE 


PHOTON 


TABLES 


TABLE  1 


OSCILLATOR  STRENGTH  FOR  PHOTOIONIZATION  FROM  1S-LEVEL 
NE-LIKE  AL:  Is2  2£2  2p^ 

ELECTRON  DENSITY=2 . 2E+23  TEMP=50  EV 


PHOTON  ENERGY (RYD) 
0.01 
0.02 
0.03 
0.04 
0.05 
0.10 
0.20 
0.50 
1.00 
2.00 
5.00 
10.00 
12.00 
15.00 
20.00 
30.00 
50.00 
100.00 


OSC . STRENGTH ( 1/RYD ) 
0.2729E-02 
0.1053E-01 
0.1829E-01 
0.1983E-01 
0.1810E-01 
0.1084E-01 
0.7484E-02 
0.5958E-02 
0.5477E-02 
0.5295E-02 
0.5059E-02 
0.4698E-02 
0.4563E-02 
0.4371E-02 
0.4069E-02 
0.3530E-02 
0.2682E-02 
0.1485E-02 


TABLE  2 


OSCILLATOR  STRENGTH  FOR  PHOTOIONIZATION  FROM  2S-LEVEL 
NE-LIKE  AL:  1sZ  2  S2*  2p^ 

ELECTRON  DENSITY=2 . 2E+23  TEMP=50  EV 


PHOTON  ENERGY (RYD) 
0.01 
0.02 
0.03 
0.04 
0.05 
0.10 
0.20 
0.50 
1.00 
2.00 
5.00 
10.00 
12.00 
15.00 
20.00 
30.00 
50.00 
100.00 


OSC . STRENGTH ( 1/RYD ) 
0.9451E-02 
0.3644E-01 
0.6331E-01 
0.6857E-01 
0.6257E-01 
0.3740E-01 
0.2569E-01 
0.2018E-01 
0.1814E-01 
0.1679E-01 
0.1419E-01 
0.1098E-01 
0.9983E-02 
0.8713E-02 
0.7046E-02 
0.4814E-02 
0.2552E-02 
0.8175E-03 


TABLE  3 


OSCILLATOR  STRENGTH  FOR  PHOTOIONIZATION  FROM  2P-LEVEL 

O.  (3 

NE-LIKE  AL:  IS  2S  2p 
ELECTRON  DENSITY=2 . 2E+23  TEMP=50  EV 


PHOTON  ENERGY (RYD) 
0.01 
0.02 
0.03 
0.04 
0.05 
0.10 
0.20 
0.50 
1.00 
2.00 
5.00 
10.00 
12.00 
15.00 
20.00 
30.00 
50.00 
100.00 


OSC . STRENGTH ( 1/RYD ) 
0.1027E-02 
0.1101E-02 
0.1185E-02 
0.1366E-02 
0.1683E-02 
0.6747E-02 
0.3191E-01 
0.2787E-01 
0 . 2268E-01 
0.2009E-01 
0.1563E-01 
0.1068E-01 
0.9262E-02 
0.7537E-02 
0.5464E-02 
0. 3092E-02 
0.1237E-02 
0.2478E-03 


TABLE  4 


OSCILLATOR  STRENGTH  FOR  PHOTOIONIZATION  FROM  1S-LEVEL 

2.  2  5 

AL  ION  :  Is  2s  2p 
ELECTRON  DENSITY=2 . 2E+23  TEMP=50  EV 


PHOTON  ENERGY ( RYD ) 
0.01 
0.02 
0.03 
0.04 
0.05 
0.10 
0 . 20 
0.50 
1.00 
2.00 
5.00 
10 . 00 
12.00 
15.00 
20.00 
30.00 
50.00 
100.00  . 


OSC . STRENGTH ( 1/RYD ) 
0 . 1006E-02 
0 . 3188E-01 
0.5733E-02 
0 .7858E-02 
0 . 9296E-02 
0 . 1008E-01 
0.7977E-02 
0.6076E-02 
0 . 5545E-02 
0.5315E-02 
0 . 5069E-02 
0.4704E-02 
0 . 4569E-02 
0 . 4376E-02 
0.4075E-02 
0.3538E-02 
0 . 2689E-02 
0 . 1 487E-02 


TABLE  5 


OSCILLATOR  STRENGTH  FOR  PHOTOIONIZATION  FROM  2S-LEVEL 
AL  ION  :  Is2"  2S2  2.p 
ELECTRON  DENSITY=2 . 2E+2 3  TEMP=50  EV 


PHOTON  ENERGY (RYD) 
0.01 
0.02 
0.03 
0.04 
0.05 
0.10 
0.20 
0.50 
1.00 
2.00 
5.00 
10.00 
12.00 
15.00 
20.00 
30.00 
50.00 
100.00  . 


OSC . STRENGTH ( 1/RYD ) 
0 . 3482E-02 
0.1103E-01 
0.1982E-01 
0.2716E-01 
0 . 3211E-01 
0.3473E-01 
0.2737E-01 
0 . 2056E-01 
0.1835E-01 
0 . 1684E-01 
0 . 1420E-01 
0 . 1 098E-01 
0.9976E-02 
0 . 8704E-02 
0 .7038E-02 
0.4809E-02 
0 . 2  550E-02 
0 . 8165E-03 


TABLE  6 


OSCILLATOR  STRENGTH  FOR  PHOTOIONI ZATION  FROM  2P-LEVEL 
AL  ION  :  /S'1  2S2  2p5 

ELECTRON  DENSITY=2 . 2E+2 3  TEMP=50  EV 


PHOTON  ENERGY (RYD) 
0.01 
0.02 
0.03 
0.04 
0.05 
0.10 
0.20 
0.50 
1.00 
2.00 
5.00 
10.00 
12.00 
15.00 
20.00 
30.00 
50.00 
100.00  . 


OSC . STRENGTH ( 1/RYD ) 
0 . 1787E-02 
0 . 1538E-02 
0.1462E-02 
0 . 1504E-02 
0.1652E-02 
0 . 4203E-02 
0 . 1878E-01 
0 . 2823E-01 
0 .2303E-01 
0 . 2015E-01 
0 . 1564E-01 
0 . 1070E-01 
0 . 9277E-02 
0 . 7551E-02 
0 . 5476E-02 
0.3097E-02 
0 . 1237E-02 
0 . 2477E-03 


TABLE  7 


OSCILLATOR  STRENGTH  FOR  PHOTOIONIZATION  FROM  1S-LEVEL 
AL  ION:  |s22S^2-p4' 

ELECTRON  DENSITY=2 . 2E+23  TEMP=50  EV 

PHOTON  ENERGY(RYD)  OSC. STRENGTH(1/RYD) 


0.02 

0.7035E-02 

0.03 

0.7789E-02 

0.04 

0.8168E-02 

0.05 

0. 8333E-02 

0.10 

0.8253E-02 

0-20 

0. 7312E-02 

0.50 

0.5549E-02 

1.00 

0. 5595E-02 

2.00 

0. 5368E-02 

5.00 

0.5075E-02 

10.00 

0.4704E-02 

12.00 

0.4569E-02 

15.00 

0.4376E-02 

20.00 

0.4075E-02 

30.00 

0. 3540E-02 

50.00 

0. 2694E-02 

100.00 

0. 1488E-02 

TABLE  8 


OSCILLATOR  STRENGTH  FOR  PHOTOIONIZATION  FROM  2S-LEVEL 
AL  ION:  \ S1  2S2  2p* 

ELECTRON  DENSITY=2 . 2E+23  TEMP=50  EV 


PHOTON  ENERGY(RYD) 

OSC. STRENGTH (1/RYD) 

0.02 

0.2409E-01 

0.03 

0.2665E-01 

0.04 

0. 2794E-01 

0.05 

0. 2849E-01 

0.10 

0. 2815E-01 

0.20 

0. 2483E-01 

0.50 

0. 1858E-01 

1.00 

0. 1833E-01 

2.00 

0. 1683E-01 

5.00 

0. 1407E-01 

10.00 

0. 1086E-01 

12.00 

0.9877E-02 

15.00 

0.8617E-02 

20.00 

0.6967E-02 

30.00 

0.4767E-02 

50.00 

0. 2531E-02 

100.00 

0.8114E-03 

TABLE  9 


OSCILLATOR  STRENGTH  FOR  PHOTOIONIZATION  FROM  2P-LEVEL 
AL  ION: 

ELECTRON  DENSITY=2 . 2E+23  TEMP=50  EV 


PHOTON  ENERGY(RYD) 

OSC. STRENGTH(1/RYD) 

O.Ol 

0. 6903E-03 

0.02 

0. 1471E-02 

0.03 

0.3150E-02 

0.04 

0.5769E-02 

0.05 

0.8491E-02 

0.10 

...  0. 1799E-01 

0.20 

0.2966E-01 

0.50 

0. 2515E-01 

1.00 

0. 2264E-01 

2.00 

0.2001E-01 

5.00 

0. 1544E-01 

10.00 

0. 1057E-01 

12.00 

0.9169E-02 

15.00 

0.7469E-02 

20.00 

0.5422E-02 

30.00 

0. 3069E-02 

50.00 

0. 1225E-02 

100.00 

0. 2461E-03 

TABLE  10 


OSCILLATOR  STRENGTH  FOR  PHOTOIONIZATION  FROM  1S-LEVEL 

AL  ION  :  IS2  2.S2  2p6  3S 

ELECTION  DENSITY=2 . 2E+23  TEMP=50  EV 


PHOTON  ENERGY ( RYD ) 
0.01 
0.02 
0.03 
0.04 
0.05 
0.10 
0.20 
0.50 
1.00 
2.00 
5.00 
10.00 
12.00 
15.00 
20.00 
30.00 
50.00  ■ 


OSC . STRENGTH ( 1 /RYD ) 
0.5854E-02 
0.6898E-02 
0.7021E-02 
0 . 6924E-02 
-•  0 . 6786E-02 
0 . 6211E-02 
0.5818E-02 
0 . 5641E-02 
0.5402E-01 
0.5278E-01 
0.5050E-01 
0.4688E-01 
0.4551E-02 
0.4355E-02 
0 . 4049E-02 
0 . 3  505E-02 
0.2666E-02 


100.00 


0.1478E-02 


TABLE  11 


OSCILLATOR  STRENGTH  FOR  PHOTOIONIZATION  FROM  2S-LEVEL 
AL  ION  :  IS2  IS**  2-p^  3S 
ELECTRON  DENSITY=2 . 2E+23  TEMP=50  EV 


PHOTON  ENERGY ( RYD ) 
0.01 
0.02 
0.03 
0.04 
0.05 
0.10 
'  0.20 
0.50 
1.00 
2.00 
5.00 
10.00 
12.00 
15.00 
20.00 
30.00 
50.00 
100.00 


OS  C . STRENGTH ( 1/RYD ) 
0 . 2009E-01 
0 . 2366E-01 
0 . 2407E-01 
0.2373E-01 
0 . 2325E-01 
0 .2123E-01 
0 . 1979E-01 
0 . 1894E-01 
0 . 1774E-01 
0.1660E-01 
0 . 1407E-01 
0 . 1091E-01 
0.99 21E-02 
0 . 8659E-02 
0 . 7003E-02 
0 . 4785E-02 
0 . 2149E-02 
0 . 8159E-03 


TABLE  12 


OSCILLATOR  STRENGTH  FOR  PHOTOIONIZATION  FROM  2P-LEVEL 
AL  ION:  IS2  2.S2  3S 
ELECTRON  DENSITY=2 . 2E+23  TEMP=50  EV 


PHOTON  ENERGY (RYD) 

OSC . STRENGTH ( 1 /RYD) 

O.Ol 

0. 2995E-03 

0.02 

0.6383E-03 

0.03 

0. 1251E-02 

0.04 

0. 2479E-02 

0.05 

0.4795E-02 

0.10 

...  0.6495E-01 

0.20 

0.4336E-01 

0.50 

0.2520E-01 

1.00 

0. 2204E-01 

2.00 

0. 1988E-01 

5.00 

0. 1543E-01 

10.00 

0. 1051E-01 

12.00 

0.9106E-02 

15.00 

0.7410E-02 

20.00 

0.5381E-02 

30.00 

0. 3060E-02 

50.00 

0. 1229E-02 

100.00 

0. 2468E-03 

TABLE  13 


OSCILLATOR  STRENGTH  FOR  FREE-FREE  TRANSITONS 
NEON-LIKE  ALUMINUM 

ELECTRON  DENSITY=2 . 2E+23  TEMP=50  EV 
PHOTON  ENERGY=0- 370  AT. UNIT 


E(INCIDENT) 

(AT. UNIT) 

E(OUTGOING) 

(AT. UNIT) 

OSCILLATOR  STRENGTH  * 
(AT. UNIT) 

1.00 

1.370 

0.1980E+01 

2.00 

2.370 

0.8505E+02 

5.00 

5.370 

0. 2601E+02 

10.00 

10.370 

0.5506E+02 

15.00 

15.370 

0. 9203E+02 

20.00 

20.370 

6. 1351E+03 

25.00 

25.370 

0.1830E+03 

30.00 

30.370 

0.2352E+03 

50.00 

50.370 

0. 2991E+03 

70.00 

70.370 

0.3707E+03 

100.00 

100.370 

0.4499E+03 

150.00 

150.370 

0.5368E+03 

200.00 

200.370 

0.6281E+03 

*  Sum  of  1>(1-1)  osc.  strengths 


TABLE  14 


OSCILLATOR  STRENGTH  FOR  FREE-FREE  TRANSITONS 
NEON-LIKE  ALUMINUM 


ELECTRON  DENSITY=2 . 2E+23 

PHOTON  ENERGY=0. 370 

TEMP=50  EV 

AT. UNIT 

E( INCIDENT) 

E( OUTGOING) 

OSCILLATOR  STRENGTH 

(AT. UNIT) 

(AT. UNIT) 

(AT. UNIT) 

1.00 

1.370 

0. 1935E+02 

2.00 

2.370 

0.4743E+02 

5.00 

5.370 

0. 9232E+02 

10.00 

10.370 

0. 1555E+03 

15.00 

15.370 

0.2325E+03 

20.00 

20.370 

0.3208E+03 

25.00 

25.370 

0.4189E+03 

30.00 

30.370 

0.5254E+03 

50.00 

50.370 

0.6574E+03 

70.00 

70.370 

0.8069E+03 

100.00 

100.370 

0.9744E+03 

150.00 

150.370 

0. 1160E+04 

200.00 

200.370 

0. 1357E+04 

Sum  of  1>(1+1)  osc.  strengths 


TABLE  15 


OSCILLATOR  STRENGTH  FOR  FREE-FREE  TRANSITONS 
NEON-LIKE  ALUMINUM 

ELECTRON  DENSITY=2 . 2E+23  TEMP=50  EV 
PHOTON  ENERGY=0. 740  AT. UNIT 


E( INCIDENT) 

(AT. UNIT) 

E(OUTGOING) 

(AT. UNIT) 

OSCILLATOR  STRENGTH  * 
(AT. UNIT) 

1.00 

1.740 

0. 1917E+00 

2.00 

2.740 

0. 7512E+00 

5-00 

5.740 

0. 2469E+01 

10.00 

10.740 

0.5598E+01 

15.00 

15.740 

0.9739E+01 

20.00 

20.740 

0. 1467E+02 

25.00 

25.740 

0. 2024E+02 

30.00 

30.740 

0.2635E+02 

50.00 

50.740 

0.3402E+02 

70.00 

70.740 

0.4269E+02 

100.00 

100.740 

0.5237E+02 

150.00 

150.740 

0. 6304E+02 

200.00 

200.740 

0.7429E+02 

*  Sum  of  1>(1-1)  osc.  strengths 


TABLE  16 


OSCILLATOR  STRENGTH  FOR  FREE-FREE  TRANSITONS 
NEON-LIKE  ALUMINUM 

ELECTRON  DENSITY=2 • 2E+23  TEMP=50  EV 
PHOTON  ENERGY=0 .740  AT. UNIT 


E(INCIDENT) 

(AT. UNIT) 

E( OUTGOING) 

(AT. UNIT) 

OSCILLATOR  STRENGTH 
(AT. UNIT) 

1.00 

1.740 

0.4742E+01 

2.00 

2.740 

0. 1021E+02 

5.00 

5.740 

0. 1734E+02 

10.00 

10.740 

0. 2646E+02 

15.00 

15.740 

0.3713E+02 

20.00 

20.740 

0.4911E+02 

25-00 

25.740 

0. 6223E+02 

30.00 

30.740 

0.7634E+02 

50.00 

50.740 

-  0.9347E+02 

70.00 

70.740 

0. 1127E+03 

100.00 

100.740 

0.1340E+03 

150.00 

150.740 

0. 1576E+03 

200.00 

200.740 

0. 1825E+03 

*  Sum  of  1>(1+1)  osc.  strengths 


TABLE  17 


OSCILLATOR  STRENGTH  FOR  FREE-FREE  TRANSITONS 
NEON-LIKE  ALUMINUM 

ELECTRON  DENSITY=2 . 2E+23  TEMP=50  EV 

PHOTON  ENERGY=1 . 295  AT. UNIT 


E( INCIDENT) 

(AT. UNIT) 

E( OUTGOING) 

(AT. UNIT) 

OSCILLATOR  STRENGTH  * 
(AT. UNIT) 

1.00 

2.295 

0. 4104E-01 

2.00 

3.295 

0. 1223E+00 

5.00 

6.295 

0. 3603E+00 

10-00 

11.295 

0.8347E+00 

15.00 

16.295 

0. 1495E+01 

20.00 

21.295 

0: 2304E+01 

t 

25.00 

26.295 

0-3238E+01 

30.00 

31,295 

0.4277E+01 

50-00 

51.295 

~  0.5621E+01 

70.00 

71.295 

0.7164E+01 

100.00 

101.295 

0.8907E+01 

150.00 

151.295 

0.1085E+02 

200.00 

201.295 

0. 1291E+02 

*  Sum  of  1>(1-1) 


osc . 


strengths 


TABLE  18 


OSCILLATOR  STRENGTH  FOR  FREE-FREE  TRANSITONS 
NEON-LIKE  ALUMINUM 

ELECTRON  DENSITY=2 . 2E+23  TEMP=50  EV 
PHOTON  ENERGY=1 . 295  AT. UNIT 


E(INCIDENT) 

E(OUTGOING) 

OSCILLATOR  STRENGTH  * 

(AT. UNIT) 

(AT. UNIT) 

(AT. UNIT) 

1.00 

2.295 

0. 1552E+01 

2.00 

3.295 

0. 3140E+01 

5.00 

6.295 

0.4921E+01 

10.00 

11.295 

0. 6979E+01 

15.00 

16.295 

0.9278E+01 

20.00 

21.295 

0. 1179E+02 

25.00 

26.295 

O'.  1448E+02 

30.00 

31.295 

0.1734E+02 

50.00 

51.295 

0.2072E+02 

70.00 

71.295 

0.2446E+02 

100.00 

101.295 

0.2857E+02 

150.00 

151.295 

0.3306E+02 

200.00 

201.295 

0.3778E+02 

*  1>(1+1)  osc.  strengths 


TABLE  19 


OSCILLATOR  STRENGTH  FOR  FREE-FREE  TRANSITONS 
NEON-LIKE  ALUMINUM 

ELECTRON  DENSITY=2. 2E+23  TEMP=50  EV 
PHOTON  ENERGY=2. 220  AT. UNIT 


E(INCIDENT) 

E( OUTGOING) 

OSCILLATOR  STRENGTH 

(AT. UNIT) 

(AT. UNIT) 

(AT. UNIT) 

1.00 

3.220 

0.4791E+00 

2.00 

4.220 

0. 9546E+00 

5.00 

7.220 

0. 1449E+01 

10.00 

12.220 

0. 1976E+01 

15.00 

17.220 

0. 2535E+01 

20.00 

22.220 

0.3126E+01 

25.00 

27.220 

0. 3745E+01 

30.00 

32.220 

0.4391E+01 

50.00 

52.220 

0.5123E+01 

70.00 

72.220 

'  0.5915E+01 

100.00 

102.220 

0. 6772E+01 

150.00 

152.220 

0. 7696E+01 

200.00 

202.220 

0.8659E+01 

*  1>(1+1)  osc.  strengths 


TABLE  20 


OSCILLATOR  STRENGTH  FOR  FREE-FREE  TRANSITONS 
NEON-LIKE  ALUMINUM 

ELECTRON  DENSITY=2 . 2E+23  TEMP=50  EV 
PHOTON  ENERGY=2 . 220  AT. UNIT 


E( INCIDENT) 

(AT. UNIT) 

E( OUTGOING) 

(AT. UNIT) 

OSCILLATOR  STRENGTH 
(AT. UNIT) 

1.00 

3.220 

0. 1293E-01 

2.00 

4.220 

0. 2996E-01 

5.00 

7.220 

0. 6580E-01 

10.00 

12.220 

0. 1374E+00 

15.00 

17.220 

0. 2416E+00 

20.00 

22.220 

0. 3740E+00 

25.00 

27.220 

0.5308E+00 

30.00 

32.220 

0. 7088E+00 

50.00 

52.220 

0.9500E+00 

70.00 

72.220 

0. 1233E+01 

100.00 

102.220 

0. 1559E+01 

150.00 

152.220 

0 . 1929E+01 

200.00 

202.220 

0. 2324E+01 

*  Sum  of  1>(1-1)  osc.  strengths 
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SAMPLE  COMPUTER  OUTPUT 


STEP  STEP1 

0. 310000E-02  0. 500000E-03  0.130000E+02 


USING  THE  DFT/INPUT  POTENTIAL 


O.OOOOOOE+OO 
0. 310000E-02 
0. 620000E-02 
0. 930000E-02 
0.124000E-01 
0. 155000E-01 
0. 186000E-01 
0. 217000E-01 
0. 248000E-01 
0. 279000E-01 
0 . 310000E-01 
0. 341000E-01 
0. 372000E-01 
0 . 403000E-01 
0. 434000E-01 
0. 465000E-01 
0.496000E-01 
0. 527000E-01 
0 . 558000E-01 
0. 589000E-01 


0. 130001E+02 
0. 129925E+02 
0. 129839E+02 
0. 129743E+02 
0. 129636E+02 
0. 129521E+02 
0. 129398E+02 
0. 129267E+02 
0. 129128E+02 
0. 128983E+02 
0. 128831E+02 
0 . 128673E+02 
0. 128510E+02 
0. 128341E+02 
0. 128168E+02 
0. 127990E+02 
0. 127808E+02 
0. 127622E+02 
0. 127433E+02 
0- 127240E+02 


Ne-like  Al  ion  in  Al  plasma 
n  2.2E+23  T=  50  ev 
Bound-Free  cross-sections 
from  lS-level 


N=  1  L=  0 


ALPHA  = 

-0. 77000000E+02 

LAST  ZERO 

POINT  I  =  261 

ZR  = 

0.28 

ALPHA  = 

-0. 77099998E+02 

LAST  ZERO 

POINT  I  =  261 

ZR  = 

0.28 

ALPHA  = 

-0. 77199997E+02 

LAST  ZERO 

POINT  I  =  261 

ZR  = 

0.28 

ALPHA  = 

-0. 77299995E+02 

LAST  ZERO 

POINT  I  =  271 

ZR  = 

0.33 

ALPHA  = 

-0. 77399994E+02 

LAST  ZERO 

POINT  I  -  271 

ZR  = 

0.33 

ALPHA  = 

-0. 77499992E+02 

LAST  ZERO 

POINT  I  =  271 

ZR  = 

0.33 

ALPHA  = 

-0. 77599991E+02 

LAST  ZERO 

POINT  I  =  271 

ZR  = 

0.33 

ALPHA  = 

-0. 77699989E+02 

LAST  ZERO 

POINT  I  =  271 

ZR  = 

0.33 

ALPHA  = 

-0. 77799988E+02 

LAST  ZERO 

POINT  I  =  271 

ZR  = 

0.33 

ALPHA  =  -0 . 77899986E+02 

LAST  ZERO  POINT  I  =  271 


ZR 


0.33 


ALPHA  = 

-0. 77999985E+02 

■  LAST  ZERO 

POINT  I  =  271 

ZR  = 

0.33 

■  ALPHA  = 

-0 . 78099983E+02 

m  LAST  ZERO 

POINT  I  =  271 

ZR  = 

0.33 

ALPHA  = 

-0. 78199982E+02 

ffl  LAST  ZERO 

POINT  I  =  271 

ZR  = 

0.33 

1  ALPHA  = 

-0. 78299980E+02 

LAST  ZERO 

POINT  I  =  271 

ZR  = 

0.33 

m  ALPHA  = 

-0. 78399979E+02 

■  LAST  ZERO 

POINT  I  =  271 

ZR  = 

0.33 

"  ALPHA  = 

-0. 78499977E+02 

LAST  ZERO 

POINT  I  =  271 

ZR  = 

0.33 

■  ALPHA  = 

-0. 78599976E+02 

■  LAST  ZERO 

POINT  I  =  281 

ZR  = 

0.38 

ALPHA  = 

-0. 78699974E+02 

■m  LAST  ZERO 

POINT  I  =  281 

ZR  = 

0.38 

9  ALPHA  = 

-0. 7879997 3E+02 

■  LAST  ZERO 

POINT  I  =  281 

ZR  = 

0.38 

ALPHA  = 

-0. 78899971E+02 

9  LAST  ZERO 

POINT  I  =  281 

ZR  = 

0.38 

9  ALPHA  - 

-0. 78999969E+02 

LAST  ZERO 

POINT  I  =  281 

ZR  = 

0.38 

_  ALPHA  = 

-0. 79099968E+02 

9  LAST  ZERO 

POINT  I  =  291 

ZR  = 

0.43 

■  ALPHA  = 

-0. 79199966E+02 

LAST  ZERO 

POINT  I  =  301 

ZR  = 

0.48 

■  ALPHA  = 

-0. 79299965E+02 

|  ALPHA  = 

-0. 79249969E+02 

LAST  ZERO 

POINT  I  =  321 

ZR  = 

0.58 

—  ALPHA  = 

-0. 79274971E+02 

g  ALPHA  = 

-0. 79262474E+02 

■  ALPHA  = 

-0. 79256226E+02 

LAST  ZERO 

POINT  I  =  341 

ZR  = 

0.68 

■  ALPHA  = 

-0. 79259354E+02 

|  ALPHA  = 

-0. 79257790E+02 

ALPHA  = 

-0. 79257011E+02 

_  ALPHA  = 

-0. 79256622E+02 

H  LAST  ZERO 

POINT  I  =  351 

ZR  = 

0.73 

■  ALPHA  = 

-0. 79256821E+02 

LAST  ZERO 

POINT  I  =  361 

ZR  = 

0.78 

■  ALPHA  = 

-0. 79256920E+02 

9  ALPHA  = 

-0. 79256874E+02 

MATCH 1  =  251 

MATCH2  = 

267 

_  ALPHA  =  -0. 79256874E+02  RATIO  -  1.0  =  -0. 180399E-02 

■  NUCLEAR  CHARGE  =  13.0  Z  =  1.0  DIPOLE  POLARIZIBILITY  =  0.0000 

I  SEMI-EMPIRICAL  RADIAL  WAVE  FUNCTION 

N  =  1  L  =  0  ENERGY  (ATOM.  UNITS)  =  -0. 7925687E+02 


DEBYE  LENGTH 


0. 100000E+01 


ATDIM  =  0.000000E+00 


0.00000 

0.00000E+00 

0.03500 

0 . 20933E+01 

0.22500 

0. 11238E+0 

0.00050 

0. 46877E-01 

0.03900 

0. 22139E+01 

0.25000 

0.90107E+0 

0.00100 

0. 93146E-01 

0.04300 

0. 23167E+01 

0.27500 

0. 71534E+0 

0.00150 

0. 13881E+00 

0.04700 

0. 24034E+01 

0.30000 

0. 56325E+0 

0.00200 

0. 18389E+00 

0.05100 

0 . 24751E+01 

0.32500 

0. 44045E+0 

0.00250 

0. 22837E+00 

0.05500 

0. 25333E+01 

0.35000 

0. 34240E+0 

0.00300 

0. 27226E+00 

0.05900 

0 . 25791E+01 

0.37500 

0. 26483E+0 

0.00350 

0. 31558E+00 

0.06300 

0. 26137E+01 

0.40000 

0. 20394E+0 

0.00400 

0. 35833E+00 

0.06700 

0 . 26380E+01 

0.42500 

0. 15644E+0 

0.00450 

0 . 40051E+00 

0.07100 

0 . 26531E+01 

0.45000 

0. 11959E+0 

0.00500 

0 . 44212E+00 

0.07500 

0. 26597E+01 

0.50000 

0. 69270E-0 

0.00600 

0. 52369E+00 

0.08000 

0 . 26574E+01 

0.55000 

0. 39729E-0 

0.00700 

0 . 60307E+00 

0.08500 

0 . 26448E+01 

0.60000 

0. 22603E-0 

0.00800 

0 . 68031E+00 

0.09000 

0 . 26231E+01 

0.65000 

0. 12774E-0 

0.00900 

0. 75545E+00 

0.09500 

0 . 25935E+01 

0.70000 

0. 71793E-0 

0.01000 

0.82853E+00 

0.10000 

0. 25572E+01 

0.80000 

0. 22370E-0 

0.01100 

0.89960E+00 

0.10500 

0 . 25151E+01 

0.90000 

0. 68658E-0 

0.01200 

0. 96869E+00 

0.11000 

0. 24680E+01 

1.00000 

0 . 20818E-0 

0.01300 

0. 10358E+01 

0.11500 

0. 24169E+01 

1 . 10000 

0 . 62524E-0 

0.01400 

0. 11011E+01 

0. 12000 

0 . 23623E+01 

1 . 20000 

0. 18637E-0 

0.01500 

0 . 11645E+01 

0.12500 

0 . 23050E+01 

1.40000 

0. 16268E-0 

0.01700 

0. 12858E+01 

0.13500 

0. 21842E+01 

1 . 60000 

0. 13944E-0 

0.01900 

0. 14001E+01 

0.14500 

0 . 20585E+01 

1 . 80000 

0. 11790E-0 

0.02100 

0 . 15076E+01 

0.15500 

0 . 19309E+01 

2.00000 

0. 98670E-0 

0.02300 

0 . 16087E+01 

0. 16500 

0 . 18037E+01 

2 . 20000 

0.81929E-1 

0.02500 

0 . 17036E+01 

0.17500 

0 . 16787E+01 

2 . 40000 

0. 67618E-1 

0.02700 

0.17925E+01 

0. 18500 

0. 15574E+01 

2.60000 

0. 55542E-1 

0.02900 

0. 18757E+01 

0.19500 

0. 14406E+01 

2.80000 

0. 45448E-1 

0.03100 

0. 19534E+01 

0.20500 

0 . 13292E+01 

3.00000 

0. 37075E-1 

0.03300 

0. 20259E+01 

0.21500 

0. 12235E+01 

3 . 20000 

0. 30067E-1 

EIGVAL 

ENERGY 

INTEGRAL 

XSECT  (PF/pe) 

-0. 792569E+02 

0. 100000E+03 

0. 735737E-02 

0. 148477E-02 

-0 . 792569E+02 

0. 500000E+02 

0. 110105E-01 

0. 268212E-02 

-0. 792569E+02 

0.300000E+02 

0. 1 32855E-01 

0.353044E-02 

-0 . 792569E+02 

0 . 200000E+02 

0. 146585E-01 

0. 406986E-02 

-0. 792569E+02 

0. 150000E+02 

0. 1 54090E-01 

0. 437135E-02 

-0. 792569E+02 

0. 120000E+02 

0. 158819E-01 

0 . 456345E-02 

-0. 792569E+02 

0. 100000E+02 

0. 162093E-01 

0. 469780E-02 

-0. 792569E+02 

0. 500000E+01 

0. 170761E-01 

0. 505896E-02 

-0. 792569E+02 
-0. 792569E+02 
-0. 792569E+02 
-0. 792569E+02 
-0. 792569E+02 
-0. 792569E+02 
-0. 792569E+02 
-0. 792569E+02 
-0. 792569E+02 
-0. 792569E+02 


0.200000E+01 
0. 100000E+01 
0. 500000E+00 
0. 200000E+00 
0. 100000E+00 
0. 500000E-01 
0.400000E-01 
0. 300000E-01 
0. 200000E-01 
0. 100000E-01 


0 . 176328E-01 
0. 179889E-01 
0 . 187920E-01 
0. 210821E-01 
0. 253847E-01 
0. 328018E-01 
0 . 343316E-01 
0. 329822E-01 
0. 250204E-01 
0 . 127382E-01 


0.529522E-02 
0. 547693E-02 
0.595815E-02 
0. 748465E-02 
0. 108446E-01 
0. 181021E-01 
0. 198288E-01 
0. 182995E-01 
0. 105303E-01 
0. 272924E-02 


LISTING  OF  SAMPLE  COMPUTER  COPES: 
FFNEW.FOR;  DIFF.FOR;  SIMP.FOR;  WAVEF.FOR. 


o  o  o  o  o 


HCGRAM  FFNEIW 


Calculates  photoabsorpticn  cross  sections  for  ocntinuun  states  in 
the  DFT  potential 

CCWN  U(1751)  ,FU(1751)  ,FV(1751)  ,FKL(10, 1751)  ,AFKL(10) , 

#  BFKL(10), 

1  Y(1751), 01(851), 0(851), P(3, 851), ENt(2),Z, SUP, ABC, EN(3), 

2  LFKL(  10) ,  XCHEOC(  10) 

OmN/QC2/  LP0T,ICKr,AMFL 
(MW  /SCREEN/  CRnX851),^PID(851) 

HMEIKTCN  ENRG(6)  ,R(840)  ,RVR(840) 

C 

C  Sun  of  1>1+1  nHtrix  elsrent  is  ffosc 

C  FF1N  contains  input  info  for  ocntiruun  states  etc 
C  KJEB  contains  the  EFT  potaitial 

C  NLNERICAL  IITOVAlEIVE  of  potaitial  is  used  in  rratrix  eLanait 
C  FFDLTT  contains  prcgran  output  of  f-f  cross  section 

C  FFOSC  is  osc.strargth  in  at.init 

C 
C 

c  0PEN(lIOT4XB,FIlB='FTlN.D!i^^,  ,SIMU5='0LD'  ) 

c  oPN(um,4]09,Fiifi='FCJiB.cm''  jSFxnjfc'aD' ) 

c  om(um'^,niE=,Fmn:.iw\simjs=,mj') 

c 

BEAD  (5,13)  ICHEEK 
13  FOMT  (1011) 

READ(5,40)SIEP,SIEP1,2; 

WRIIE(06,4001) 

4001  KFMAT(8X,  'SUP'  ,8X,  'STEP1' ) 

WRTIE(06,4000)SIEP,SIEP1,ZZ 
4CXX)  F0BMAT(1X,3EB.6) 

40  KFmr(3E33.6) 

VKHE(6,2050) 

2050  FCFMAT(//, '  USIFG  HE  CFT/IMUr  FOENIJAL' ,//) 

A=6450.0*SEEP1 

Z=1.0 

IQKT=0 

U(1)=0.0 

LF0r=O 

IfcSEP 

J=1 

Fcao 

108  DO  109  I=J,K 
U(I+1)4J(J)+(I-J+1)*H 

109  ONIINLE 

IF  (K-250)  110,111,2000 

110  H=2*SIEP 
J=51 
K=250 
GOTO  108 

111  H=10*SIEP 
J=251 
K=1750 
GOTO  108 

2000  0CNEIRE 


10  FCFMAT  (1E13.6) 

C 

C  VREAD  reads  in  r  and  rV(r)  vectors 
C 

CALL  VREAD(R,roR) 

C 

C  fill  GRID  and  YC3HD  for  interpolation 
C 

DO  20  1=1,840 
GRID(I)=£(I) 

20  7GRID(I)=KVR(I) 

DO  21  1=841,851 
GRID(I)=0.0 

21  1GRID(I)=0.0 
DO  22  1=1,851 

22  Q(I)=U(I) 

CALL  WAVEF(1.0,SIEP,CR,Q) 

DO  23  1=1,851 
GRID(I)=U(I) 

23  Y(I)=CR(I) 

DO  2  1=851,1751 
2  Y(I)=0.0 
C 

C  Ein  is  the  incident  electron  eiergy.Ebut  the  gritted  energy  in  at  .unit 
C 

SlML=0.0 

17  READ  (5,10)EIN 
CMECA  =14.80 
B0UT=EIN  +  CMECA 
DO  14  LL=0,35 
LLP1  =  ILfl 
CALL  DUF(LL,E3N,1) 

CALL  nnT(UPl,EDUT,2) 

CR(1)=0. 

CR(2)=0. 

CR(851)=0. 

DO  100  1=3,850 

YEA  =  Y(I)/U(I)-Y(I-1)/U(I-1) 

YEB  =  Y(I+1)/U(I+1)-Y(I-1)/U(I-1) 

XEA  =  U(I)-U(I-1) 

XEB  =  U(I+1)-U(I~1) 

XI2  =  (YDA*>nj-YII^>LA)/(>LA*}LB*(>IlV-X[E) ) 

XII  =  YLA/XCA  -  XL2*XDA 
CR(I)  =  XU  +  2.*XEA*XE2 
100  CR(I)=ai(I)*FKL(l,I)*FKL<2,I) 

CALL  S3MP(X) 

A=XW2/CMEGA**3 
XLLPIdIPl 
XLL  =  LL 

FTDGC=(16.OtXLLPl*A)/3.0 

C 

C  Sun  over  all  1  values 
C 

SIML  =SLML*-FPD9C 
14  CENTIME 

WRTIE  (6,51) 

51  FCFMAT(//,5X,  'CMECA'  ,5X,  'E3N'  ,10X,  'EMIT' ,  1CK,  '+SJML' ) 


o  o  o  o  o  o  o 


52 


WRITE  (6,52)CMEEA,EIN,HXJr,Sl«L 
KFMfflT(/,4H3.6) 

©TO  17 
END 


SLEROUTOE  VREAD(R,KvR) 
nD®EKNR(840),  FMR(840) 

READS  EFT  EFTEUUVE  RHENHAL  IMl 

READ(9,11)(R(I),I=1,840) 

11  KFMAlT(6E13.6) 

READ(9,ll)(RVR(I),I=l,&W) 
WRnE(06,33)((R(I),RyR(T)),I=l,100,4) 
33  FDBMA[r(lX,E13.6,2X,EI3.6) 

WRTIE(6,106) 

106  RHWT(///) 

REMN 


EM) 


SlHOjmE  DEF(L,ENC,IFKL) 


c 

C  Numerical  Solution  Of  Ccntinjun  States  of  Schrodinger  Bq. 

C 

ocmw  11(1751)  ,FU(1751)  ,FV(1751)  ,FKL(10, 1751)  ,AFKL(10)  ,BFKL(10) , 

1  Y(1751),CR(851), 0(851), P(3, 851), EN1(2),Z, SUP, ABC, EN(3), 

2  LFKL(10),iafrK(10) 
amwcc2/  iPor,icRr,AMPL 
HMEM3CN  FAC(150) 

FEALA8  BL,ENER,li\,lB,UC,AU,BU,CU,EU,Bl,B2,B3,B4,B5,B6 

KEAL*8  BUS, GE, DUG 

LFKL(IFKL)=L 

IF  (L.IT.0)  FEMN 

EFER=EMC*l.QDtOO 

AL=L 

IF  (EMR.IT. 0.0)  GO  TO  1325 

IF  (IIOT.EjQ.0.AM).ENER.IT.lD-03)  GD  TO  6010 

EMR=FMF/Z**2 

ST=1.0 

J1=0 

00=0.0 

01=0.0 

IF  (L.IT.10)  GO  TO  6006 
IF  (ICKT.NE.O)  WtTIE  (6,6013)  L 

6013  RFMAT  (////, 5X.14HWAFNIIG  L  =  ,I3,/5X,37HFKL  FCR  L.GT.9  ARE  NOT 
1  CB3HDQGNAIIZED,  //// ) 

IGKT=0 

Y1=Y(17)-Y(1) 

Y2=Y(33)-Y(1) 

Z1=(4*YL-Y2)/(32*SEEP) 

Z2=(Y2-2*YL)/(512*SIEP**2) 

K=51 

GOTO  6007 

6006  K.il 
Y1=Y(6)-Y(1) 

Y2=Y(11)-Y(1) 

Z2=(4*Y1-Y2)/(10SIEP) 

Z2=(Y2-2*Y1)/(50*SIEP**2) 

IF  (IOFT.EQ.O)  GO  TO  6007 
A=0(6)/U(6)**(Lfl) 

BdQ(ll)/U(ll)**(Lt-l) 

QO=2*A-B 

IF  (QO.IT.O.O)  00=0.0 
Q1=(B-Q0)/U(11) 

6007  EB=>-Y(l)*Sr/(ALfl.O) 

CG=(  -2*Y(l)*EB-2*Zl*Sr-2*EMR*Sr)/(4*AL+6.0) 

GOGCTO0/(4*AL*6.0) 

CD=(  -2*Y(l)^0C-2*ZL*EB-2*Z2*ST-2*Einy®)/(6*'ALfl2 .0) 

ED=O>Ql/(6*ALfl2.0) 

FU(1)=0.0 

B=1.0 

1246  DO  1245  1=2, K 
A=U(I)*B 

IF  (L.ED.0)  GO  TO  1245 


AsA*U(I) 

IF  (ABS(A) .  LI.  IE-20)  A=0.0 
6001  OCNTIME 

1245  FU(I)=(ST+BB^(I)+aXJ(I)^+0>U(I)**3)*A 
BL=L*(Lfl) 

IF  (ABS(FU(K)).IT.lE-20)  0)  TO  6008 
IF  (20-K)  302,302,6009 

6008  J1=J1+1 
B=10B 

GO  TO  1246 

6009  H=SIEP**2 
K=50 
Kl=10 
K2=ll 
K3=12 

300  Bltl.0-(BL/U(K1)**2-2.0D>Y(K1)/IJ(K1)-2.0*®ER''*H/I2.0 
atl.0-(HyU{K2)^-2.OX>Y(K2)/U(K2)-2.0^ER)*H/12.0 
EU=l.O-(BLOJ(IG)^-2.aX^Y(K3)AJ(K3)-2.0*®FR)W12TO 
ALtl2.0-10.0*CU 
UA=HJ(KL)*1.0D0 
IB=FU(K2)*1 .  CDO 
BUG=UDOQ(Kl) 
a&=inxKK2) 

0JG=UXXXK3) 

DO  301  I=K2,K 
lXt(ALKB-JBU*UA)/DLI 
IF  (ICKT.BQ.O  )  GO  TO  6012 
IF  (I.GT.850)  GO  TO  6012 
imi>(BLDfl(XlJ0+aE)*{l/12.0 
6012  FU(I+1)=UC 

IF  (DABS(UC) .  IT.  lEf25)  GO  TO  3050 

J1=J1-10 

IB=LB*1D-10 

UC4JOTO-10 

JM=I+1 

5000  FU(JM)=FU(JM)*1E-10 
IF  (JM.00.2)  GO  TO  3050 
IF  (ABS(FU( JM) ) .  LI.  IE-20)  GO  TO  5002 
JK=JM-1 
GOTO  5000 

5002  DO  5003  KFW2.JM 

5003  FU(KM)=0.0 
3050  0LK3J 

emu 

DU=1 .0-(BL/U(I+2)**2-2  .CD0*Y(I+2)/U(I+2)-2 . 0*®IR)*H/12 .0 

Alfcl2.0-10.0*CU 

UA=DB 

IF  (ICKT.H3.0.0R.I.GT.849)  GO  TO  301 

EUKJUG 

aG=ai; 

amiaxxi+2) 

301 IB=UC 

IF  (K-250)  302,303,304 
302  H=4*SIEP**2 
10*49 
K2=51 
K3=52 


£=250 
GOTO  300 

303  fklOO*SIEP**2 
KL=246 
K2=251 
K3=252 
£=1749 
GOTO  300 

304  RJ(1751)=(ALWU(1750)-BLRI4)/EU 
J=0 

DO  3049  1=10,1005 

IF  ( (FU(I) .IT.0.0.AND.FU(I+1) .GT.O.O) .CR. (FU(I) .GT.O.O.M). 

1  FU(I+1).IT.0.0)>  J=J+1 

3049  OCNIIRE 
K=1 

MAMPW) 

DO  305  1=1006,1748 

if  «FU(i).ir.o.o.M).nj(i+i).Gr.o.o).®.(HJ(i).Gr.o.o.m 
1  RJ(I+1).LT.0.0))  GO  TO  306 
GOTO  999 

306  Bl=FU(I-2)*l.dWD0 
B2=FU(I-1)-B1 
B&=(FU(I)-Bl-2*B2)/2.0 
B4=(FU(I+l)-Bl-3*B2-6*B3)/6.0 
B5=(FU(I+2)-Bl-4^B2-12*B3-24^34)/24.0 
B6=(FU(I+3)-Bl-5*B2-20"ffl-60*B4-120*B5)/ 120.0 
B4J(I)+10»SIEP*ABS(FU(I))/(ABS(FU(I))4ABS(RJ(I+1))) 

D1=B-  U(I-2) 

D2=B-  U(I-1) 

D3=&-U(I) 

D4=&-  U(I+1) 

D6=B-  U(I+2) 

D6=B-  U(  1+3) 

G=0.1/SIEP 

A=C  ^+CK>®3*(Dl+ffi)+C^3*B4*(Dl^+Dl;M034E:TO)4>A4  *B5* 

1  (Dl*D2mmWI4^^  W(Dl*D2*D3*D4f 

2  DMMB*C6TOM2^^ 

D=Y(I)+C  *(B4J(I))*(Y(I+1)-Y(I)) 

O^I7B^*2+2.0*D/Bf2.0*EMR 

IF  (C.LT.O.O)  MAMPL=1 
IF  (C.LT.O.O)  GO  TO  999 
0=ABS(C) 

C=3CKT(C) 

C=3TT(C) 

F/C(K)=ABS(Cy(AWr(Z) ) ) 

IF  (ICHErK(l).NE.O)  WRTIE  (6,70Q)  L,K,B,F«XK) 

703  FOWiT  (2I5,F10.4,E20.8) 

K=K+1 

IF  (K.B0.71)  GO  TO  3040 
999  OONUNLE 
306  0CNI3MJE 
3042  A=G.O 
K=K-1 

IF  (MAMPL.HU)  J=0 

IF  (MAMFL.H3.1.M).K.H).0)  4RIIE  (6,6011)  L,  ®RG 
IF  (K.EQ.0)  GO  TO  1220 
DO  307  1=1, K 


307  A=A+F/C(I) 

WC 

IF  (J.ED.O)  A=F£C(K) 

AMPLfA 
GOTO  1221 

1220  WRITE  (6,1222)  L 

AMPLTO.O  _ 

1222  FORMAT  (1(K,29HAMPLTIIIE  NOT  EEEMMD  L  =,  15) 
GOTO  1325 

1221  DO  308  1=1,1750 
J=1752-I 
B=A*FU(J) 

IF  (ABS(B)  .GT.  IE-08)  GO  TO  6000 
B=FU(J-1)*A 

IF  (ABS(B) .  IT.  IE-08)  GO  TO  3046 
6000  OCMINUE 

308  FU(J)=B 
GOTO  3047 

3046  DO  3048  1=1,  J 
3048  FU(I)=0.0 

3047  aMTRE 
B=15450*SIEP 

A^(4i/B^2+2*Y(1751)/B+2*fl®)*5W 

A=ABS(A) 

A=SCFT(A) 

A=SCET(A) 

A=A*FU(1751) 

IF  (ABS(A).GT.0.9)  GO  TO  1210 
GOTO  1212 
1210  1=1751 

4001  1=1-2 

IF  (FU(I)/FU(1751).LT.0.71)  GD  TO  4000 
GOTO  4001 

4000  A=2.0*rLFOT/U(I  )-BLAI(I  )**2+2*H d 
A=3Jir(A) 

Bl=9CPr(A)*SCFT(Z) 

0=0.0 

J=I 

4002  &^.0*lJGr/U(I+l)^UU(I+l)**2+2*EFBR 
B=3CFT(B) 

0=2.0MJGT/U(I+2)-BLTO(I+2)^+2*EIEK 

C=3FT(C) 

eh>a+4*m; 

1=1+2 

A=C 

IF  (I.U.1751)  GO  TO  4002 

D=O*TO.0*SIEP/3.0 

A=B1*FU(J) 

0=1.0 

IF  (FU(J+1  LT.FU(J))  0-1.0 
&=C*SCFT'(J  .0-A*A) 

0=5IN(P) 

D=COo(D) 

aFKUIFKL^TO+BKJ 
BFKL(IIKL>W>-A*C 
GO  TO  4000 


1212  0=1.0 

IF  (FU(1751).IT.FU(1750))  C=-1.0 
B=C*9TT(1.0A*A) 

1213  AFKL(IIKL)=A 
BFKL(IFKL)=£ 

4003  00  3052  1=1,1751 
3052  FKL(1IKL,I)=FU(I) 

FETORN 

1325  DO  1326  1=1,1751 

1326  FU(I)=0.0 
A=0.0 
0=0.0 
GOTO  1213 

3040  AAX^€R*Z**2 

IF  (L.EQ.O)  WEE  (6,3041)  AAX 

3041  FCFMAT  (//,1CK,2CHENEFGY  (AT.  IMIS)  =,E15.6/1(X,97HEK)GRAM  MAY  NO 
IT  GIVE  dFRECT  AMPUME  OF  aMTOULM  EINZICN5  AT  SUCH  HIGH  ENRG 
2Y  -  DETFEASE  SIEP,//) 

GOTO  3042 

6010  VRTIE  (6,6011)  L,EMC 

6011  KFMAT  (/,  10K,52WAFMN0  -  HttRAM  OMDT  CCMPUIE  EKL  RR  3m  EN 
1FFGY/1CK,4HL  =  ,I5,lCK,Sffi«CY  =  ,0.4.6) 

GOTO  1325 


SIHWJmE  UAVEF(Z,STEP,Q,U) 
nMJGICN  0(851), U(851),FU(1751),FV(1751) 

D3MEMIICN  NP(14)  ,ZET(14)  ,CLEM(14)  ,IFL(1751) 

cnttN  /smw  amx85i),iraaD(85i) 

C  KMAX  MAY  EE  SET  BJM,  TO  14  IF  MFEER  OF  TIFMS  IN  THE  EXPANSION 
C  FOR  THE  RADIAL  WAVE  FINUCN  IS  GEEAER  THAN  7 
KMAX=7 

READ  (5,1)  M,(NP(I), 1=1, KMAX) 

IF  (M.GT.O)  GO  TO  2 

1  KRMAT  (1512) 

READ  (5,3)  (ZET(I), 1=1, KMAX) 

READ  (5,3)  (OFM(I), 1=1, KMAX) 

3  FORMAT  (7F10.5) 

DO  23  1=1, KMAX 
te£*NP(I) 

A=2*ZET(I) 

B=A 

IF  (M.BQ.O)  GO  TO  24 
DO  25  MF=1,M 

25  A=A/MF 

CLFM(I)=OFM(I)*B«NP(I) 

24  (UEM(I)=a£M(I)*SCFr(A) 

23  OCNUNLE 
Q(1)=0.0 
DO  26  1=2,851 
R=U(I)/Z 
0=0.0 

DO  27  J=1,KMAX 

IF  (NP(J).EQ.O)  GO  TO  31 

A=R*2ET(J) 

IF  (A.GT.23.0)  GO  TO  28 
WXP(-4) 

GOTO  29 

28  BD.O 

29  G=C^BARHNP(J)*C1EM(J) 

27  OUnMJE 

31  Q(I)=C 

26  (3M3NJE 
REHJFN 

2  IF  (M.ED.2)  GO  TO  101 
IF  (M.FD.l)  GO  TO  112 
DO  113  1=1,851 

113  FV(I)=7CRID(I) 

DO  115  1=852,1751 
115  FV(I)=TORID(851) 

112  OCNimE 
J=1 
K=10 

PMAX=6450*SEEP/Z 

11  IF  (M.ED.l)  READ  (5,5)  (FU(I),I=J,K) 

IF  (M.FD.l)  GO  TO  110 
DO  111  I=J,K 

111  FiKi)=amxi) 

110OM3ME 
5  FOFMAT  (10F8.5) 

IF  (FU(K).IT.FMAX)  GO  TO  6 


WBTIE  (6,7) 

7  FCFMAX  (//,KK,62MAR®G  -  WAVE  FUNCTION  EMEDS  AUUMD  FH3CN  - 
1  MPEAS  SEP) 

STOP 

6  IF  (K.LE.1740)  GO  TO  8 
wrte:  (6,9) 

9  FORMAT  (//,1CK,68E1WAFNIIC  -  WAVE  FUNCTICN  GIVEN  BY  MFBER  OF  VALLE 
IS  QREAER  THAN  1740) 

STOP 

8  IF  (FU(K).GT.O.O)  03  TO  10 
INHC-10 

GOTO  12 

10  J=J+10 
FGK+10 
GOTO  11 

12  DO  13  1=1,  IN 

13  FU(I)=FU(I)*Z 
A=FU(IN)-FU(IN-1) 

DO  14  I=IN,1750 

14  FU(I+1)=FU(I)4A 

IF  (FU(1751)  .GT.  6450*SIEP)  G3  TO  103 
FU(1751)  =  6450*SIEP  +  5.0 
103  DO  100  I  =  1,  1751 

100  TPL(I)  =  FU(I) 

GO  TO  4 

101  DO  102  I  =  1,  1751 

102  FU(I)  =  TPL(I) 

4  IF  (M.BQ.3)  GO  TO  114 
READ  (5,5)  (FV(I),I=1,IN) 

DO  15  I  =  IN,  1750 

15  FV(I+1)  =  0.0 
114  CONTINUE 

0(1)  =  0.0 

DO  16  I  =  2,  851 

DO  17  J  =  1,  1751 

IF  (FU(J)  .GT.  U(I))  GO  TO  18 

17  OCNEINLE 

18  A  =  FU(J+1)  -  FU(J  -  1) 

B  =  FU(J)  -  FU(J  -  1) 

C  =  FU(J+1)  -  FU(J) 

Q(I)=FV(J-1)+((FV(J)-FV(J-1))*'A*(FU(J+1)-U(I))+(FV(J+1)-FV(J-1)) 

1  *B*(U(I)-FU(J)))*(U(I)-FU(J-1))/ (A*B*C) 

16  ouiinue 

IF  (M.FQ.3)  Q(1)=2*Q(2)-Q(3) 

RETUFN 

ETC 


SUBROUTINE  SIMP (THETA) 

COMMON  U( 1751 ), FU( 1751 ),  EV( 1751 ),FKL( 10, 1751 )  ,AFKL( 10 ),BFKL( 10), 

1  Y( 1751 ) ,QR(851) ,Q(851 ) ,P( 3,851) ,ENR(2) ,Z,STEP,ABC,EN(3) , 

2  LFKL( 10 ) , ICHECK( 10 ) 

H=STEP 

THETA=H* ( 4 . 0*QR{ 2 ) +QR( 3 )  )/3 . 0 

JA=3 

KA=49 

801  DO  802  IA=JA, KA, 2 

802  THETA=THETA+H*(QR(IA)+4*QR(IA+l)+QR(IA+2)  )/3.0 
IF  ( KA-249 )  803,804,805 

803  H=2*STEP 
JA=51 
KA=249 

GO  TO  801 

804  H=10*STEP 
JA=251 
KA=849 

GO  TO  801 

805  RETURN 
END 


